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ABSTRACT 



While the origin of r-process nuclei remains a long-standing mystery, recent 
spectroscopic studies of extremely metal-poor stars in the Galactic halo strongly 
suggest that it is associated with core-collapse supernovae. In this study we 
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examine r-process nucleosynthesis in a "prompt supernova explosion" from an 
8 — IOMq progenitor star, as an alternative scenario to the "neutrino wind" mech- 
anism, which has also been considered to be a promising site of the r-process. 
In the present model, the progenitor star has formed an oxygen-neon-magnesium 
core (of mass 1.38Mq) at its center. Its smaller gravitational potential, as well 
as the smaller core that is in nuclear statistical equilibrium at the time of core 
bounce, as compared to the iron cores in more massive stars, may allow the 
star to explode hydrodynamically, rather than by delayed neutrino heating. The 
core-collapse simulations are performed with a one-dimension, Newtonian hydro- 
dynamic code. We obtain a very weak prompt explosion, in which no r-processing 
occurs. We further simulate energetic prompt explosions by enhancement of the 
shock-heating energy, in order to investigate conditions necessary for the produc- 
tion of r-process nuclei in such events. The r-process nucleosynthesis is calculated 
using a nuclear reaction network code including relevant neutron-rich isotopes 
with reactions among them. The highly neutronized ejecta {Yg ^ 0.14 — 0.20) 
leads to robust production of r-process nuclei; their relative abundances are in 
excellent agreement with the solar r-process pattern. Our results suggest that 
prompt explosions of 8 — lOM© stars with oxygen-neon- magnesium cores can be a 
promising site of r-process nuclei. The mass of the r-process material per event is 
about two orders of magnitude larger than that expected from Galactic chemical 
evolution studies. We propose, therefore, that only a small fraction of r-process 
material is ejected, owing to the "mixing-fallback" mechanism of the core matter, 
wherein most of the r-process material falls back onto the proto-neutron star. 

A lower hmit on the age of the universe is derived by apphcation of the U- 
Th chronometer pair by comparison with the observed ratio of these species in 
the highly r-process enhanced, extremely metal-poor star CS 31082-001. The 
inferred age is 14.1 ±2.4 Gyr - the same as that obtained previously based on 
the neutrino wind scenario with the same nuclear mass formula. This suggests 
that chronometric estimates obtained using the U-Th pair are independent of the 
astrophysical conditions considered. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: abun- 
dances — supernovae: general 
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1. INTRODUCTION 

The astrophysical origin of the rapid neutron-capture (r-process) species has been a long- 
standing mystery. Recently, however, a number of important new clues have been provided 
by spectroscopic studies of extremely-metal-poor stars in the Galaxy. The appearance of 
neutron-capture elements in these oldest stars in the Galaxy, including the pure-r-process 
origin of elements such as thorium and uranium, strongly suggests that the r-process nuclei 
have come from core-collapse supernovae (Sneden et al. 1996, 2000, 2003; Gayrel et al. 
2001; Gowan et al. 2002: Hill et al. 2002). Ishimaru & Wanajo (1999) have shown that 
the large star-to-star dispersion of the observed abundances of neutron-capture elements 
relative to iron in very metal-poor stars is also naturally explained if the r-process elements 
originate from a limited mass range of core-collapse supernovae with little iron production 
(8 - IOM0 or > 3OM0). Qian & Wasserburg (2003; see also Wasserburg & Qian 2000; Qian 
& Wasserburg 2001, 2002) have proposed that the production of the heavy r-process nuclei 
{A > 130) is decoupled from the production of iron-peak and a nuclei by comparison of 
observed abundances among the extremely metal-poor stars. In this view, the production 
site of the heavy r-process nuclei is associated with the accretion-induced collapse (AIG) of 
a (carbon-oxygen or oxygen-neon-magnesium) white dwarf in a binary system (Nomoto & 
Kondo 1991), or Type II supernovae from 8 — IOMq stars (Nomoto 1984). 

So far, the "neutrino wind" scenario, in which the free nucleons accelerated by the intense 
neutrino flux near the neutrino sphere of a core-collapse supernova assemble to heavier nuclei, 
has been believed to be the most promising astrophysical site of the r-process (Woosley et 
al. 1994). Even this scenario, however, encounters some difficulties (Qian & Woosley 1996; 
Hoffman et al. 1997; Cardall & Fuller 1997; Otsuki et al. 2000; Thompson, Burrows, & 
Meyer 2001; Wanajo et al. 2001, 2002). For example, Wanajo et al. (2001) have shown that 
an extremely compact proto-neutron star, e.g., 2.OM0 and 10 km, must be formed in order 
to account for the solar r-process pattern, at least within the framework of a spherically 
symmetric explosion. Although its possibility cannot be excluded, such a compact remnant 
is allowed by only a few of the many existing equations of state of nuclear matter. 

In addition, recent spectroscopic studies of extremely metal-poor stars in the Galactic 
halo indicate that the observed abundance patterns of the hghter (Z < 56) and heavier 
{Z > 56) neutron-capture elements cannot be explained by a single astrophysical site (e.g., 
neutrino winds); there must exist at least two different r-process sites (Ishimaru Sz Wanajo 
2000; Qian & Wasserburg 2001; Johnson & Bolte 2002; Sneden & Gowan 2003). Hence, it is 
of special importance to consider alternative possibilities for the occurrence of the r-process 
in core-collapse supernovae. 

Nomoto (1984, 1987) has shown that 8— IOMq stars form an electron-degenerate oxygen- 
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neon-magnesium (0-Ne-Mg) core that does not undergo further nuclear burning; rather, it 
directly collapses due to electron capture on ^^Mg and ^°Ne (Miyaji et al. 1980). Hillebrandt, 
Nomoto, & Wolff (1984) have demonstrated that the collapsing 0-Ne-Mg core explodes in 
a prompt manner, and Wheeler, Cowan, & Hillebrandt (1998) have suggested that the 
exploding 0-Ne-Mg core could be a viable site for the r-process. It has been pointed out 
that, if the core exploded hydrodynamically prior to the onset of delayed neutrino heating 
(i.e., it underwent a prompt explosion), the electron fraction (electron number per baryon). 
Ye, in the innermost layer of the ejecta would approach ~ 0.2 (Hillebrandt, Nomoto, & Wolff 
1984). Earlier works have in fact shown that a robust r-process proceeds in such conditions 
(Schramm 1973; Sato 1974; Hillebrandt, Takahashi, & Kodama 1976). 

Recently, Sumiyoshi et al. (2001) have demonstrated that the prompt explosion of an 
IIMq star with an iron core might also be a promising site of the r-process. They obtained 
a prompt explosion by an adiabatic core-collapse calculation without inclusion of electron 
capture and neutrino transport. Thompson, Burrows, & Pinto (2003) have shown, however, 
that no explosion is obtained with such a progenitor (Wooslcy & Weaver 1995), when includ- 
ing electron capture along with a detailed treatment of neutrino transport. Many previous 
works have suggested that even the lowest mass of core-collapse supernovae (~ IOA'/q), in 
which a relatively smaller iron core is formed, may have difficulties in achieving a hydrody- 
namic explosion (Bowers & Wilson 1982; Burrows & Lattimer 1983, 1985; Bruenn 1989a,b; 
Baron & Cooperstein 1990). 

On the other hand, the question of whether 8 — IOA'/q stars that form 0-Ne-Mg cores can 
explode hydrodynamically is still open. The possibility that these stars explode promptly 
remains because of the smaller iron core present at the onset of the core bounce, as well as 
the smaller gravitational potential of their collapsing cores. Hillebrandt, Nomoto, & Wolff 
(1984) have obtained a prompt explosion of a 9Mq star with a I.SSMq 0-Ne-Mg core, while 
others, using the same progenitor, have not (Burrows & Lattimer 1985; Baron, Cooperstein, 
& Kahana 1987). Mayle & Wilson (1988) obtained an explosion, not by a prompt shock, but 
by late-time neutrino heating. Similar results can be seen in the studies of AlCs. Note that 
an AlC is an analogous phenomenon to a collapsing 0-Ne-Mg core resulted from a single 
8 — IOMq star, since both consist of electron degenerate cores (Nomoto & Kondo 1991). 
Fryer et al. (1999) have obtained an explosion by neutrino heating, while others have not 
(Baron et al. 1987; Woosley & Baron 1992). The reason for these different outcomes is due, 
perhaps, to the application of different equations of state for dense matter, although other 
physical inputs may also have some influence (Pryer et al. 1999). Thus, even if a star of 
8 — IOMq exploded, it would be difficult to derive, with confidence, the physical properties 
as well as the mass of the ejected matter. Given this highly uncertain situation it is necessary 
to examine the resulting r-process nucleosynthesis in explosions obtained with different sets 
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of input physics. 

The purpose of this study is to investigate conditions necessary for the production of r- 
process nuclei obtained in purely hydrodynamical models of prompt explosions of collapsing 
0-Ne-Mg cores, and to explore some of the consequences if those conditions are met. The core 
collapse and the subsequent core bounce are simulated by a one-dimensional hydrodynamic 
code with Newtonian gravity (§ 2). For simplicity, neutrino transport is not taken into 
account. Hence, we focus on only hydrodynamical explosions just after core bounce, not on 
the delayed explosions obtained via late-time neutrino heating. As seen in § 2, the explosion 
is marginal, and no r-processing is expected. In order to obtain r-processed material, we find 
it necessary to force the occurrence of more energetic explosions. The energetic explosions are 
simulated by artificial enhancements of the shock-heating energy, rather than by application 
of different sets of input physics, for simplicity. The r-process nucleosynthesis in these 
explosions is then calculated with the use of a nuclear reaction network code (§ 3). The 
resulting contribution of the r-process material created in these simulations to the early 
chemical evolution of the Galaxy is discussed in § 4. The results of chronometric age dating, 
using the U-Th chronometer pair based on our nucleosynthesis results, is discussed in § 5. 
A summary and conclusions follow in § 6. 



2. PROMPT EXPLOSION 

A pre-supcrnova model of a 9Mq star is taken from Nomoto (1984), which forms a 
1.38 Mq 0-Ne-Mg core near the end of its evolution (see also Miyaji et al. 1980; Miyaji & 
Nomoto 1987; Nomoto 1987). We link this core to a one- dimensional implicit Lagrangian 
hydrodynamic code with Newtonian gravity (Bowers & Wilson 1991). This core is modeled 
with a finely zoned mesh of 200 mass shells (2 x IO'^Mq to 0.8Mq, 5 x lO'^M© to 1.3Mq, 
and 5 x 10~^ - 1 x to the edge of the core). 

The equation of state of nuclear matter (EOS) is taken from Shen et al. (1998), which is 
based on relativistic mean field theory. The equation of state for the electron and positron 
gas includes arbitrary relativistic pairs as well as arbitrary degeneracy. Electron and positron 
capture on nuclei, as well as on free nucleons, are included, along with the use of the up- 
to-date rates from Langanke & Martinez-Pinedo (2000). The capture is suppressed above 
the neutrino trapping density, taken to be 3 x 10^^ g cm~^, since the neutrino transport 
process is not taken into account in this study. This simplification may be appropriate, since 
neutrinos at the early epoch of the core bounce do not appear to make a significant correction 
to the entropy compared to that obtained from shock heating as shown by Hillebrandt, 
Nomoto, & Wolff (1984). The delayed neutrino heating may not significantly modify the 
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mass trajectories of the outgoing matter either, since the bulk of the ejecta are hfted to 
~ 1000 km during the first few 100 ms from the onset of the core bounce (as shown below), 
which is far from the location of the neutrino sphere (a few tens of kilometers). It is evident, 
however, that an accurate treatment of neutrino transport will be needed to obtain accurate 
mass trajectories in future work. 

Nuclear burning is implemented in a simplified manner. The composition of the O- 
Ne-Mg core is held fixed until the temperature in each zone reaches the onset of oxygen- 
burning, taken to be 2 x 10^ K (Nomoto 1987), at which point the matter is assumed to be 
instantaneously in nuclear statistical equilibrium (NSE). The temperature is then calculated 
by including its nuclear energy release. It should be noted that we find a weak a- rich 
freezeout in the subsequent post-processing nucleosynthesis calculations (§ 3), owing to the 
entropy, ~ lOA^^/c, in the ejecta. This means that the outgoing ejecta are not in perfect 
NSE, which is assumed in our hydrodynamic calculations. An improvement that takes the 
non-NSE matter properly into account will be needed to obtain more accurate trajectories. 

We begin the hydrodynamical computations with this pre-supernova model, which has 
a density of 4.4 x 10^*^ g cm~^ and temperature of 1.3 x 10^*^ K at its center. The inner 
O.IM0 has already burned to NSE. As a result, the central Yg is rather low, 0.37, owing 
to electron capture. The core bounce is initiated when ~ 90 ms has passed from the start 
of the calculation. At this time the NSE core contains only I.OM0, which is significantly 
smaller than the cases of collapsing iron cores (> 1.3Mq). The central density is 2.2 x 
lO^'' g cm~^, significantly lower than that of Hillcbrandt, Nomoto, & Wolff (1984, slightly 
above 3 x 10^"^ g cm^^), although the temperature (= 2.1 x 10^° K) and Yg (= 0.34), are 
similar. This difference is perhaps due to the use of a relatively stiff EOS in this study. 

We find that a very weak explosion results, with an ejected mass of 0.008Mq and an 
explosion energy of 2 x 10^^ ergs (model QO in Table 1). The time variations of the radius, 
temperature, and density of each zone arc displayed in Figure 1. The lowest in the 
outgoing ejecta is 0.45, where no r-processing is expected given the entropy of ~ lOA^^^- 
This is in contrast to the very energetic explosions, with ejected masses of O.2M0, explosion 
energies of 2 x 10^^ ergs, and low Yg of ~ 0.2 obtained by Hillebrandt, Nomoto, & Wolff 
(1984). This might be a consequence of the lower gravitational energy release owing to the 
EOS applied in this study. 

In order to examine the possible operation of the r-process in the explosion of this 
model, we artificially obtain explosions with typical energies of ~ 10^^ ergs by application of 

a multiplicative factor (/shock) to the shock-heating term in the energy equation (models Q3, 
Q5, and Q6 in Table 1). We take this simplified approach in this study, since the main 
difference between our result and that by Hillebrandt, Nomoto, & Wolff (1984) appear to 
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be the lower central density in ours. If the inner core reached a higher density at the time 
of core bounce by applying, for example, a softer EOS, the matter would obtain higher 
shock-heating energy. This is clearly not a self-consistent approach, and a further study is 
needed to conclude whether such a progenitor star explodes or not, taking into account a 
more accurate treatment of neutrino transport, as well as with various sets of input physics 
(like EOSs). It should be emphasized, however, that our purpose in this paper is not to 
justify the prompt explosions of collapsing 0-Ne-Mg cores, but to investigate the conditions 
necessary for the production of r-process nuclei from such an event if it occurs. 

Tabic 1 lists the multiplicative factor applied to the shock-heating term (/shock), ex- 
plosion energy (Eexp), ejected mass (Mgj), and minimum Ye in the cjccta obtained for each 
model. Energetic explosions with E^xp > 10^^ ergs are obtained for /shock > 1.5 (models Q5 
and Q6), in which deeper neutronized zones are ejected by the prompt shock, as can be seen 
in Figure 2 (model Q6). This is in contrast to the weak explosions with Eexp < 10^° ergs 
(models QO and Q3), in which only the surface of the core blows off (Figure 1). Note that 
the remnant masses for models Q5 and Q6 are I.IQMq and O.94M0, respectively, which are 
significantly smaller than the typical neutron star mass of ~ IAMq. We consider it likely 
that a mass of ~ IAMq is recovered by fallback of the once-ejected matter, as discussed in 
§4. 

In Figure 3 the electron fraction in the ejecta of each model is shown as a function of 
the ejected mass point, M^j. For models QO and Q3, decreases steeply with M^j, since 
the duration of electron capturing is long, owing to the slowly expanding cjccta (Figure 1). 
For models Q5 and Q6, on the other hand, Yg decreases gradually with Mgj, owing to the 
fast expansion of the outgoing ejecta. Nevertheless, the inner regions approach very low 
Yf., 0.30 and 0.14 for models Q5 and Q6, respectively, owing to their rather high density 
(~ 10^^ g cm~^) at the time of core bounce (Figure 2). Note that, for model Q6, Yg increases 
again for Mgj > 0.3Mq. This is due to the fact that the positron capture on free neutrons 
overcomes the electron capture on free protons when the electron degeneracy becomes less 
effective in the high temperature matter. The innermost, slowly outgoing region suffers from 
this effect. This can be also seen in the results of Mayle & Wilson (1988), who obtained a 
neutrino-powered (not prompt) explosion with the same 0-Ne-Mg core. 

The minimum Ye (~ 0.40) in Mayle & Wilson (1988) is, however, significantly higher 
than ours. This is mainly due to neutrino capture on free nucleons, since the matter is driven 
by intense neutrino flux in their simulation. Ye might increase further once nucleons begin 
to assemble into a particles and heavy nuclei (the "a-effect" ; Meyer, McLaughlin, & Fuller 
1998). In the case of the prompt explosions considered here, however, these neutrino effects 
may not alter Ye significantly. The reason is that the bulk of the ejecta are lifted to ~ 1000 km 
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at the arrival of the delayed neutrinos (a few 100 ms from the onset of the core bounce), at 
which the capture timescale of neutrinos on free nucleons is no less than a few seconds (see, 
e.g., eq. (1) in Qian et al. 1997). The dynamical timescales of the outgoing mass shells in 
model Q6 at Tg = 1,3,5, and 7 (after the core bounce), defined by Tdyn = \p/ {dp/ dt)\^ are 
shown in Figure 4. As can be seen, the dynamical timescale prior to the r-process phase 
(Tg > 3) is significantly smaller than that of neutrino interaction. 

The trend of the Y^. — Mcj relation up to Mcj ~ O.2M0 is similar in models Q5 and 
Q6, although it is inverted at Mpj ~ 0.14Mq, owing to the slightly different contribution of 
the positron and electron capture on free nucleons (Figure 3). Hence, the — Mcj relation 
between the surface and the innermost layer of the ejecta is expected to be similar to that of 
model Q6, as long as the explosion is sufficiently energetic (> 10^^ ergs). In the subsequent 
sections, therefore, we focus only on model Q6, which is taken to be representative of cases 
where r-process nucleosynthesis occurs. The ejected mass, Mej, is thus taken to be a free 
parameter, instead of simulating many other models by changing /shock- Note that the results 
by Hillebrandt, Nomoto, & Wolff (1984) and Sumiyoshi et al. (2001) are very similar to the 
cases with Mej ~ 0.2 — O.SM© in model Q6. 

3. THE r-PROCESS 

The yields of r-process nucleosynthesis species, adopting the model described in § 2 for 
the physical conditions, arc obtained by application of an extensive nuclear reaction network 
code. The network consists of ~ 3600 species, all the way from single neutrons and protons up 
to the fermium isotopes {Z = 100). Wc include all relevant reactions, i.e., (n, 7), {p, 7), {a, 7), 
{p,n), {a,p), {a,n), and their inverses. Reaction rates are taken from Thielemann (1995, 
private communication) for nuclei with Z < 46 and from Cowan, Thielemann, & Truran 
(1991) for those with Z > 47. The latter used the mass formula by Hilf et al. (1976). The 
three-body reaction a{an,^yBe, which is of special importance as the bottleneck reaction 
to heavier nuclei, is taken from the recent experimental data of Utsunomiya et al. (2001). 
The weak interactions, such as /3-decay, /3-delayed neutron emission (up to three neutrons), 
and electron capture are also included, although the latter is found to be unimportant. 

The a-decay chains and spontaneous fission processes are taken into account only after 
the freezeout of all other reactions, as in Wanajo et al. (2002). For the latter, all nuclei with 
A > 256 are assumed to decay by spontaneous fission only. The few known nuclei under- 
going spontaneous fission for A < 256 are also included, along with their branching ratios. 
Neutron- induced and /3-delayed fissions, as well as the contribution of fission fragments to 
the lighter nuclei, are neglected. Obviously, these treatments of the fission reactions are 
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oversimplified. Nevertheless, this may be acceptable, at least to first order. We leave more 
accurate treatment of these matters to future work. 

Each calculation is started at Tg = 9 (where Tg = K). The initial composition 

is taken to be that of NSE with the density and electron fraction at Tg — 9, and consists 
mostly of free nucleons and alpha particles. 

For reasons outlined in § 2, we examine the nucleosynthesis in model Q6 only, in 
which robust r-processing is expected. The resulted abundances in several representative 
Lagrangian mass shells are depicted in Figure 5. The initial electron fractions are 0.29, 0.20, 
0.18, 0.16, 0.14, and 0.20 for the corresponding zone numbers 83, 92, 95, 98, 105, and 132, 
respectively. As can be seen, a robust r-processing is possible only for Yg ^ 0.20. In partic- 
ular, a substantial amount of thorium and uranium are produced only when Yg is less than 
0.18. 

The mass-integrated abundances from the surface (zone 1) to the zones 83, 92, 95, 
98, 105, and 132 are compared with the solar r-process abundances (Kappeler et al. 1989) 
in Figure 6 (models Q6a-f in Table 2). The latter is scaled to match the height of the 
first {A = 80) and third {A = 195) peaks of the abundances in models Q6a-b and Q6c-f, 
respectively. The ejecta masses of these models are listed in Table 2. The nucleosynthesis 
result in model Q5 (not presented here) is expected to be similar to that of model Q6a, 
because of the resemblance of the Yg — -^ej profiles between these models (Figure 3). As 
can be seen in Figure 6, a solar r-process pattern for A > 130 is naturally reproduced in 
models Q6c-f, while models Q6a-b fail to reproduce the third abundance peak. This implies 
that the region with Ye < 0.20 must be ejected to account for production of the third r- 
process peak. Furthermore, to account for the solar level of thorium {A = 232) and uranium 
{A = 235,238) production, the region with rather low Yg (< 0.18) must be ejected. Note 
that large deficiencies of nuclei at A ^ 115 may be supplied if the slower rates of neutron 
capture in this region arc adopted, as demonstrated by Goriely (1997). 

We find that, for models Q6c-e, the lighter r-process nuclei with A < 130 are somewhat 
deficient compared to the solar r-process pattern (Figure 6c-e). This trend can be also 
seen in the observational abundance patterns of the highly r-process-enhanced, extremely 
metal- poor stars CS 22892-052 (Sneden et al. 2003) and CS 31082-001 (Hill et al. 2002). 
In model Q6f, the deficiency is outstanding because of large ejection of the low Ye matter 
(Figure 3). This is in contrast to the previous results obtained for the neutrino wind scenario, 
which significantly overproduce the nuclei with A 90 (Woosley et al. 1994; Wanajo et al. 
2001). The nuclei with A < 130 can be supplied by slightly less energetic explosions, like 
models Q6a-b (Figures 5a-b). It is also possible to consider that these lighter r-process nuclei 
originate from "neutrino winds" in more massive supernovae (> lOM©). The nuclei with 
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A < 130 can be produced naturally in neutrino winds with a reasonable compactness of the 
proto-neutron star, e.g., IAMq and 10 km (Wanajo et al. 2001). 

Figure 6 implies that the production of thorium and uranium differs from model to 
model, even though the abundance pattern seems to be universal between the second and 
third r-process peaks, as seen in models Q6c-f. This is in agreement with recent obser- 
vational results suggesting that the ratio Th/Eu may exhibit a star-to-star scatter, while 
the abundance pattern between the second and third peaks is in good agreement with the 
solar r-process pattern (Honda et al., in preparation). Thus, the use of Th/Eu as a cos- 
mochronometer should be regarded with caution, at least until the possible variations can 
be better quantified; U/Th might be a far more reliable chronometer, as discussed further 
in § 5. 

It is interesting to note that the nucleosynthesis results obtained for models Q6c-e are 
in good agreement with that of the prompt explosion of a IIMq star with an iron core 
by Sumiyoshi et al. (2001). Our results seem, however, in better agreement with the solar 
r-process pattern, in particular near the rare-earth peak (^4 160) and the third peak 
{A 195). This is a consequence of the slower ejection of the innermost region in our results 
owing to the reduction of pressure by electron capture, which is not taken into account in 
Sumiyoshi et al. (2001). In particular, the dynamical timescale of outgoing matter becomes 
significantly long during the epoch of the r-process (Tg ~ 3 to 1), as can be seen in Figure 4. 
The extremely low Yg (< 0.2) drives matter near the neutron-drip line, where the abundance 
pattern deviates from the solar one. In our models, however, the solar r-process pattern 
is recovered by the "freezeout effect" after the epoch of r-processing, as discussed in detail 
by Wanajo et al. (2002). This arises because the matter stays at higher temperatures for a 
longer time, hence the quasi-equilibrium between neutron emission by photodisintegration 
and subsequent neutron-capture processes continues to operate (Surman & Engel 2001). 

It is also interesting to note that the mass of the r-processed ejecta in our results 
(> O.O5M0) is more than one order of magnitude larger than that obtained by Sumiyoshi 
et al. (2001, ~ 0.003Mq). This is also a consequence of the longer dynamical timescale 
of outgoing ejecta in our simulation. The a particles are mostly consumed prior to the r- 
process phase in our calculations, whereas they dominate the final nucleosynthesis yields in 
Sumiyoshi et al. (2001). These results strongly suggest that an accurate treatment of electron 
capture, as well as of neutrino transport, will be crucial to derive an accurate prediction of 
the r-process abundances from prompt explosions of 0-Ne-Mg cores. 
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4. CONTRIBUTION TO CHEMICAL EVOLUTION OF THE GALAXY 

We have demonstrated in the previous section that, if stars of 8 — IOM0 do explode ener- 
getically, the solar r-process pattern can be reproduced quite naturally. It is of importance, 
however, to see if the nucleosynthetic material contributed to the Galaxy from such stars is 
consistent with the currently observed elemental abundances in the solar neighborhood. 

Figure 7 shows the "production factors" per supernova for model Q6e, defined for each 
nuclide as the final mass fraction, X^j, diluted by the total ejected mass (T.QM©) from the 
9Mq star, divided by its solar abundance Xq (Anders & Grevesse 1989). The solid fines 
connect isotopes of a given element (after nuclear decay). The dotted horizontal lines indicate 
a "normalization band" between the largest production factor (^^^Xe) and that by a factor 
of ten less than that, along with a median value (dashed fine). This band is taken to be 
representative of the uncertainty in the nuclear data for very neutron-rich nuclei (Woosley 
et al. 1994). As can be seen in Figure 7, a majority of the nuclei with A > 100 fall within 
the normafization band (except for the large deficiencies near A = 120), which is regarded 
to be the dominant species produced. The deficiencies near A = 88, 138, and 208 may be 
supplied by s-process contributions from other sources. Although such a supernova might be 
considered a fikely production site of r-process nuclei, its contribution to nuclei with A < 80, 
including the a elements and iron-peak species, is negligible. 

One of the essential questions raised by previous works is that, if prompt supernova 
explosions are one of the major sites of r-process nuclei, would in fact the r-process nuclei be 
significantly o?;erproo?uceo? (Hillebrandt, Takahashi, & Kodama 1976). As far as the explosion 
is purely hydrodynamical, a highly neutronized deeper region must be ejected in order for a 
successful r-process to result. It seems inevitable, therefore, that one must avoid an ejection 
of large amounts of r-process matter, at least when assuming spherical symmetry. Our 
result shows that more than O.OSM© of the r-process matter {A > 120) is ejected per event, 
which reproduces the solar r-process pattern (models Q6c-f in Table 2). This is about three 
orders of magnitude larger than the 5.8 x 10~^Mq in the neutrino-heated supernova ejecta 
from a 2OM0 star obtained by Woosley et al. (1994). This value might would be reduced 
to some extent by a-rich freezeout in faster outgoing mass shells than ours as observed in 
Sumiyoshi et al. (2001). However, good agreement with the solar r-process pattern would 
not be achieved in the ejecta expanding too fast, as discussed in § 3. 

As discussed by Woosley et al. (1994), production factors must be on the order of ~ 10 
for the case that all supernovae contribute equally to r-process production in the Galaxy. 
Stars of 8 — IOMq would account for ~ 30% of all supernovae (if they explode). The ejected 
masses (~ 7 — 9Mq) are smaller than a factor of two to three than those from more massive 
supernovae. Thus, production factors of ~ 100 are allowed in the case that all 8 — IOM0 
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stars contribute equally to the Galactic r-process material. The production factors in our 
results are, however, about three orders of magnitude higher than this (Figure 7). 

It might be argued that this type of event is extremely rare, accounting for only 0.01 — 
0.1% of all core-collapse supernovae. However, observations of extremely metal-poor stars 
([Fe/H] ~ -3) in the Galactic halo show that at least two, CS 22892-052 and CS 31082-001, 
out of about a hundred studied at high resolution, imply contributions from highly r-process- 
enhanced supernova ejecta (Sncdcn ct al. 2000; Hill ct al. 2002). Moreover, such an extremely 
rare event would result in a much larger dispersion of r-process elements relative to iron 
than is currently observed amongst extremely metal-poor stars. Ishimaru & Wanajo (1999) 
demonstrated that the observed star-to-star dispersion of [Eu/Fe] over a range ~ —1 to 2 dex, 
was reproduced by their chemical evolution model if Eu originated from stars of 8 — IOMq. 
Recent abundance measurements of Eu in a few extremely metal-poor stars with [Fe/H] < 
—3 by SUBARU/HDS further supports their result (Ishimaru et al., in preparation). The 
requisite mass of Eu in their model is ~ 1O~^M0 per event. The ejected mass of Eu in our 
result is more than two orders of magnitude larger (Table 2). 

In order to resolve this conflict, we propose that the "mixing- fallback" mechanism op- 
erates in this kind of supernova. The peculiar abundance patterns of some extremely iron- 
deficient stars, including HE 0107-5240 with [Fe/H] = -5.3 (Christheb et al. 2002), is ex- 
plained successfully with this mechanism, as proposed by Umeda & Nomoto (2002, 2003). 
If a substantial amount of the hydrogen and helium envelope remains at the onset of the 
explosion, the outgoing ejecta may undergo large-scale mixing by Rayleigh- Taylor instabili- 
ties, which is believed to have happened in SN1987A (Hachisu et al. 1990; Herant & Woosley 
1994; Kifonidis et al. 2003). Thus a tiny amount, say, ~ 1%, of the r-process material is 
mixed into the outer layers and then ejected, but most of the core material may fall back 
onto the proto-neutron star via the reverse shock arising from the hydrogen-helium layer 
interface. In this case, the typical mass of the proto-neutron star (~ I.4M0) is recovered. 
An asymmetric explosion mechanism, such as that which might arise from rotating cores, 
may have a similar effect as the ejection of deep-interior material in a small amount (Yamada 
& Sato 1994; Fryer & Heger 2000). This "mixing-fallback" scenario must be further tested 
by detailed multidimensional-hydrodynamic studies. However, it may provide us with a new 
paradigm for the nature of supernova nucleosynthesis. 

This type of event may be characterized as a "faint" supernova, owing to the weakened 
explosion energy as well as the reduced amount of ^^Ni by "mixing-fallback" . In addition, its 
ejecta consists mostly of hydrogen and helium with a higher ratio of He/H than that in the 
solar system. This event can be easily distinguished from the core-collapse supernovae of iron 
cores resulted from stars of > lOM©, which are characterized by abundant ct-elements. It is 
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interesting to note that such an explosion bears a close resemblance to the Crab supernova, 
SN1054 (Nomoto et al. 1982). There is the possibility, therefore, that direct evidence will be 
obtained from the Crab nebula (or other similar supernova remnants) through the detection 
of r-process elements by future observations (Wallcrstein ct al. 1995). Detection of the 7-ray 
lines from decays of radioactive isotopes produced by the r-process could provide additional 
direct evidence for this scenario (Qian, Vogel, & Wasserburg 1998). 

It is of special importance to confirm, from spectroscopic studies of extremely metal- 
poor stars in the Galactic halo, that the r-process elements are not associated with the 
production of a- and iron-peak elements (Wasserburg & Qian 2000; Qian & Wasserburg 
2001, 2002). Qian & Wasserburg (2003) have suggested that the r-process enrichment in 
extremely metal-poor stars HD 115444, HD 122563, and CS 31082-001 is independent of 
the production of the elements from Na to Zn (including a- and iron-peak elements). The 
abundances of Na-Zn among these stars are mostly the same, while the level of r-process 
enhancement differs from star to star. These authors take this as evidence that the heavy 
r-process nuclei originated from AICs or Type II supernovae from 8 — IOMq stars, although 
they prefer AICs for the r-process site. Note that the nucleosynthetic outcome of an AIC 
event may be similar to the collapsing 0-Ne-Mg core resulting from a single 8 — IOMq star. 
Hence, it is possible that AICs also undergo prompt explosions. The frequency of the AIC 
events in the Galaxy is expected to be very small, ~ 10~^ yr"-*^ (Bailyn Sz Grindlay 1990). 
Assuming that the frequency of the core-collapse supernovae in the Galaxy is ~ 10~^ yr~^ 
(Cappellaro et al. 1997), the production factor per AIC event is ~ 10^. This is in reasonable 
agreement with our result of 10^, as can be seen in Figure 7. In the case of AICs, however, 
there is no chance to undergo "mixing-fallback" . owing to the absence of the outer envelope. 
Thus, a much larger dispersion of [Eu/Fe] than that observed in the extremely metal-poor 
stars seems difficult to avoid. 

We consider, therefore, that the heavy r-process nuclei originate from 8 — IOMq stars 
discussed in this study, being independent of whether they are single stars or are in binary 
systems. The lighter r-process elements {Z < 56) observed in extremely metal-poor stars 
may represent simply the interstellar medium from which these stars were formed, which 
originated from, perhaps, "neutrino winds" in supernovae from stars of > IOMq. The 
observed abundances of Eu relative to iron in extremely metal-poor stars are in fact well 
reproduced with this assumption in the chemical evolution model by Ishimaru & Wanajo 
(1999). Future spectroscopic studies of extremely metal-poor stars may be able to distinguish 
the scenario suggested in this study from AICs, for example, through observations of s- 
process abundances originating from mass transfer in binary systems. 
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5. 



U-Th COSMOCHRONOLOGY 



The recently discovered r-process-enhanced, extremely metal-poor star, CS 31082-001 
([Fe/H] = —2.9) has provided a new, potentially quite powerful cosmochronometer, uranium 
(Cayrel et al. 2001; HiU et al. 2002). Wanajo et al. (2002) have determined the age of this 
star (more precisely, the time that has passed since the production event of the r-process 
species incorporated into this star) using the U-Th chronometer pair to be 14.1 ± 2.5 Gyr, 
based on the neutrino wind scenario. This is regarded as a hard lower limit of the age of 
the universe, and is in good agreement with that derived by an site-independent approach 
(Goriely & Arnould 2001; Schatz et al. 2002). It is of interest to examine if the same age of 
this star is also obtained, based on the prompt explosion model of a 9Mq star described in 
the previous sections. Note that the same mass formula (Hilf et al. 1976) as in Wanajo et 
al. (2002) is used in this study for comparison. 

In Figure 8 the available spectroscopic abundance data for CS 31082-001 (dots) are com- 
pared with the nucleosynthesis result of model Q6e discussed in § 3 (thick line) and with the 
solar r-process pattern (thin line), scaled at Eu {Z = 63). The data for the neutron-capture 
elements in this star are taken from Hill et al. (2002). An overall agreement of our result with 
the spectroscopic data up to lead {Z — 82) can be seen, although our result appears some- 
what deficient for the lighter elements. This is in contrast to the results in neutrino winds 
(Wanajo et al. 2002), in which the hghter elements are significantly overproduced. Thus, 
this model might be a reasonable one for the r-process events that produce large amounts of 
thorium and uranium, and whose products now appear in the atmosphere of CS 31082-001. 

Figure 9 shows the mass-integrated abundance ratios of Th/Eu (open squares) and 
U/Th (open circles) from the surface to the mass point Mej, as well as the inferred ages 
of CS 31082-001, t*(Th/Eu) (filled squares) and t*(U/Th) (filled circles). As discussed by 
Hill et al. (2002), the age of CS 31082-001 can be inferred by apphcation of the following 
relations: 



where the half lives of ^^^Th (14.05 Gyr) and ^^^U (4.468 Gyr), and the subscripts "0" and 
"now" denote the initial and current values derived by theory and observation, respectively. 
In principle, it*(U/Th) may serve as a more precise chronometer than t* (Th/Eu), owing to 
the smaller coefficients in front of equations (2). Moreover, the ratio U/Th is less dependent 
on the model parameter Mej, since these species are separated by only two units in atomic 
number. Note that ^^^U is assumed to have a-decayed away because of its relatively short 



r (Th/Eu) 

r(u/Th) 



46.67 [log(Th/Eu)o - log (Th/Eu) now] Gyr 
21.76 [log(U/Th)o - log(U/Th)„o^] Gyr 



(1) 
(2) 



half life (0.704 Gyr). 
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As can be seen in Figure 9, U/Th approaches a constant value (= 0.51) for Mgj > O.26M0 
{Ye < 0.17), while Th/Eu varies widely. As a result, the age of CS 31082-001 determined 
by U/Th results in a constant value, t*(U/Th) = 14.1 Gyr, for Mej > O.26M0. The age 
t*(Th/Eu) is sensitive to the parameter M^j, ranging from a negative age to 23.8 Gyr, which 
illustrates that caution must be used in the application of this chronometer pair. It is useful, 
however, to take t* (Th/Eu) as a constraint on the model parameter M^j, although the result 
might be changed if one includes an accurate neutrino transport and other input physics. 
It is found that the models with Mej = O.SOMq and 0.37Mq (zone numbers 105 and 118, 
respectively) give the same ages between T (Th/Eu) and r(U/Th). The former corresponds 
to model Q6e, which would provide a consistent scenario for the origin of the r-process 
elements in CS 31082-001, as can be seen in Figure 8. It should be noted that the fission 
fragments in model Q6e account for 23% of the mass contained in A > 120 nuclei (the last 
column in Table 2), whose contribution to the abundance pattern is neglected in this study. 
Obviously, more accurate treatment of fission reactions is needed. Nevertheless, the age 
i*(U/Th) would not be altered significantly by this improvement, since the ratio U/Th is at 
the saturated value (=0.51) at M^j — 0.26Mq, where the mass fraction of fission fragments 
is only 2%. 

It is noteworthy that the inferred age by U/Th in this study is the same as the result in 
Wanajo et al. (2002), which is based on a different astrophysical scenario (neutrino winds), 
but uses the same nuclear mass formula (Hilf et al. 1976). Our result confirms the robustness 
of the age determination using the U-Th chronometer pair, which is mostly independent of 
the astrophysical conditions considered. Rather than the r-process site, the nuclear mass 
formulae adopted, as well as the treatment of fission are crucial as far as the U-Th pair is 
concerned (Seeger & Schramm 1970; Goriely & Arnould 2001; Schatz et al. 2002). 

6. SUMMARY AND CONCLUSIONS 

We have examined the r-process nucleosynthesis obtained in the prompt explosion aris- 
ing from the collapse of a 9Mq star with an 0-Ne-Mg core. The core collapse and subsequent 
core bounce were simulated with a one-dimensional, implicit, Lagrangian hydrodynamic code 
with Newtonian gravity. Neutrino transport was neglected for simplicity. We obtained a very 
weak explosion (model QO) with an explosion energy of ~ 2 x 10"^^ ergs, and an ejected mass 
of ~ O.OO8M0. No r-processing occurred in this model, because of the high electron fraction 
(> 0.45) with low entropy (~ IOA^a^;). 

We further simulated energetic explosions by an artificial enhancement of the shock- 
heating energy, which might be expected from calculations with other sets of input physics, as 
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well as with other pre-supernova models. This resulted in an explosion energy of > 10^^ ergs 
and an ejected mass of > 0.2Mq. Highly neutronized matter {Yg 0.14) was ejected, which 
led to strong r-processing (model Q6). Material arising from r-process nucleosynthesis was 
calculated with a nuclear reaction network code containing ~ 3600 isotopes with all relevant 
reactions. The result was in good agreement with the solar r-process pattern, in particular 
for nuclei with A > 130. Some of lighter r-process nuclei {A < 130) were deficient, which is 
consistent with the abundance patterns of the highly r-process enhanced, extremely metal- 
poor stars CS 22892-052 and CS 31082-001. This implies that the fighter r-process nuclei 
may originate from another site, which we suggest might be associated with the "neutrino 
wind" in core-collapsing supernovae of iron cores resulted from more massive stars (> IOMq). 

The large ejection of r-process material (> O.OSM© per event) confiicts with the level of 
dispersion of r-process elements relative to iron observed in extremely metal-poor stars. We 
suggest, therefore, that only a small fraction (~ 1%) of the r-processed material is ejected, 
while the bulk of such material falls back onto the proto-neutron star by the "mixing-fallback" 
mechanism. 

The age of the highly r-process enhanced, extremely metal-poor star, CS 31082-001 
was derived by application of the U-Th chronometer pair, and can be regarded as a hard 
lower limit on the age of the universe. The age obtained is 14.1 ±2.4 Gyr (the quoted 
error only includes that arising from the observations), the same as that based on a difi^ercnt 
astrophysical site, the neutrino-wind scenario (Wanajo et al. 2002), using the same nuclear 
mass formula. This confirms that the age determined by the U-Th pair is mostly independent 
of the astrophysical environment considered. The dependence of age dating on different 
nuclear mass formulae, based on the prompt explosion scenario presented in this paper, will 
be reported in future work. 

It is obvious that more studies, including an accurate treatment of neutrino transport 
and other input physics, as well as multi-dimensional simulations of prompt supernova ex- 
plosions, are needed to derive the final conclusion on the r-process scenario presented in this 
study. Nevertheless, this scenario is attractive as a promising site of the r-process, since the 
solar r-process pattern can be naturally reproduced without the problematic overproduction 
of A PS! 90 that appeared in the neutrino wind scenario. This type of event is characterized 
with the absence of a- and iron-peak elements, which can be easily distinguished from that 
of the core-collapsing iron core resulted from a more massive star (> IOMq). Future spectro- 
scopic studies of extremely metal-poor stars in the Galactic halo will reveal if the collapsing 
0-Ne-Mg cores of 8 — IOMq stars are a viable site for the production of r-process nuclei. 
Improved observational determination of the U/Th ratio in CS 31082-001, as is presently 
being pursued, as well as a measured abundance of Pb in this star (as is being obtained 
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with the Hubble Space Telescope), and the identification of a greater number of highly r- 
process-enhanced, metal-poor stars, also underway, will surely deepen our understanding of 
the relevant processes involved. 
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Table 1. Results of Core-Collapse Simulations 



model 


yshock 


£;exp (10^1 ergs) 


Mej (Me) 


Y ■ 

^ e,mm 


QO 


1.0 


0.018 


0.0079 


0.45 


Q3 


1.3 


0.10 


0.029 


0.36 


Q5 


1.5 


1.2 


0.19 


0.30 


Q6 


1.6 


3.5 


0.44 


0.14 
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Table 2. Ejected Mass {Mq) 



model 




A > 120 


56 


Fe 


En 




fission"'' 


QO 


0.0079 


0.0 


0.0018 


0.0019 


0.0 




0.0 


Q6a 


0.19 


2.6 X 10-^ 


0.018 


0.020 


0.0 




0.0 


Q6b 


0.24 


0.035 


0.018 


0.020 


2.4 X 10" 


-4 


1.4 X 10-6 


Q6c 


0.25 


0.051 


0.018 


0.020 


4.1 X 10- 


-4 


4.3 X 10-^ 


Q6d 


0.27 


0.064 


0.018 


0.020 


4.3 X 10- 


-4 


0.047 


Q6c 


0.30 


0.080 


0.018 


0.020 


4.6 X 10- 


-4 


0.23 


Q6f 


0.44 


0.21 


0.018 


0.020 


0.0020 




0.15 



Wass fraction of fission fragments {A > 120) 
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Fig. 1. — Time variations of (a) radius, (b) temperature, and (c) density for all mass points in 
the weak prompt explosion of a 9Mq star (model QO). The ejected mass points are denoted 
in black, while those of the remnant are in grey. 
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Fig. 2. — Time variations of (a) radius, (b) temperature, and (c) density for selected mass 
points (with roughly an equal mass interval) in the energetic prompt explosion of a 9Mq 
star, in which the shock-heating energy is enhanced artificially by a factor of 1.6 (model Q6). 
The ejected mass points are denoted in black, while those of the remnant are in grey. Thick 
lines denote the mass shells with Mej = 0.2 and O.SMq (Ye — 0-20 and 0.14, respectively; see 
Figure 3) . The material between these lines is particularly of importance to account for the 
solar r-process abundances (see discussion in § 3). 
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Fig. 3. — Ye distribution in the ejected material in models QO (open triangles), Q3 (filled 
triangles), Q5 (open circles), and Q6 (filled circles). The surface of the 0-Ne-Mg core is at 
mass coordinate zero. For model Q6, selected mass points are denoted by zone numbers (see 
Table 2). 




Fig. 4. — The dynamical timescales of the outgoing material as functions of M^j at Tg — 
1, 3, 5, and 7 in model Q6. 
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Fig. 5. — Final abundances as a function of mass number from r-process calculations for 
trajectories (a) 83, (b) 92, (c) 95, (d) 98, (e) 105, and (f) 132 in Table 2. 
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Fig. 6. — Final mass-averaged r-process abundances (line) as a function of mass number 
obtained from the ejected zones in (a) models Q6a, (b) Q6b, (c) Q6c, (d) Q6d, (e) Q6e, and 
(f) Q6f (see Table 2). These are compared with the solar r-process abundances (points) of 
Kappclcr ct al. (1989), which is scaled to match the height of the first peak (A — 80) for 
(a)-(b) and the third peak {A = 195) for (c)-(f). 
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Fig. 7. — Mass-averaged production factors in model Q6e (see Table 2). Isotopes of a given 
element are connected by lines. Elements with even and odd atomic numbers are denoted 
by points and triangles, respectively. The dotted lines indicate a normalization band (see 
text), with its median value (dashed line). 
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Fig. 8. — Comparison of the mass-integrated yields (thick hne) for model Q6e, scaled at Eu 
{Z — 63), with the abundance pattern of CS 31082-001 (filled circles, with observational 
error bars), as a function of atomic number. For Pb, the observed upper limit is shown by 
the open circle. The scaled solar r-process pattern is shown by the thin line. 
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Fig. 9. — Mass-integrated abundance ratios Th/Eu (open squares) and U/Th (open circles) 
from the surface of the core to the mass point in model Q6. The surface of the 0-Ne- 
Mg core is at mass coordinate zero. Ages of CS 31082-001 t* (Th/Eu) (filled squares) and 
t*(U/Th) (filled circles) inferred by these ratios are also shown. The lines of these ages have 
intersections near the mass points 105 and 118. 



